Spinful Composite Fermions in a Negative Effective Field 
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In this paper, we study fractional quantum Hall composite fermion wavefunctions at filling frac- 
tions V = 2/3, 3/5, and 4/7. At each of these filling fractions, there are several possible wavefunctions 
with different spin polarizations, depending on how many spin-up or spin-down composite fermion 
Landau levels are occupied. We calculate the energy of the possible composite fermion wavefunc- 
tions and we predict transitions between ground-states of different spin polarizations as the ratio 
of Zeeman energy to Coulomb energy is varied. Previously, several experiinents have observed such 
transitions between states of differing spin polarization and we make direct comparison of our pre- 
dictions to these experiments. For more detailed comparison between theory and experiment, we 
also include finite-thickness effects in our calculations. We find reasonable qualitative agreement 
between the experiments and composite fermion theory. Finally, we consider composite fermion 
states at filling factors v = 2 + 2/3, 2 -I- 3/5, and 2 + 4/7. The latter two cases we predict to be spin 
polarized even at zero Zeeman energy. 

PACS numbers: 73.43.-f; 71.10.Pm 



I. INTRODUCTION 

The distinctive band structure of two-dimensional elec- 
trons in a magnetic field has proved to be a rich setting 
for new and exciting physical phenomena. In particu- 
lar, in the absence of interactions or disorder, the spec- 
trum of single particle eigenstates breaks into degenerate 
bands called Landau levels (LLs); fractional quantum 
Hall (FQH) physicsi occurs when interactions between 
electrons break the degeneracy of partially filled Landau 
bands, leading to an incompressible fluid ground-state. 
Quantum Hall systems are typically characterized by a 
filling factor i^, which quantifies the ratio of the number of 
electrons in the two-dimensional electron gas (2DEG) to 
the number of flux quanta through the sample or, alter- 
natively, the filling factor quantifies the number of filled 
LLs in the system. 

Naively, one might expect that the high magnetic 
field characteristic of the fractional quantum Hall regime 
might remove any spin degree of freedom entirely. How- 
ever in conventional GaAs systems, owing to the small g 
factor, even in fairly high fields, unpolarized or partially 
polarized quantum Hall states may occur. The favorabil- 
ity of spin non-polarized states is generally determined by 
the ratio of the Zeeman energy per electron, Ez = gfJ-sB 
the energy associated with flipping a spin in magnetic 
field strength B, to the Coulomb energy per electron Ec, 
the energy associated with the spatial configuration of 
electrons in the quantum Hall system. For many quan- 
tum Hall states (particularly composite fermion states) 
due to the Pauli exclusion principle, the Coulomb energy 
(or, equivalently, the effective composite fermion kinetic 
energy) can be lower if spins are not fully polarized. As 
a result, if the Zeeman energy is not too large, quantum 
Hall ground-states may not be fully spin polarized. 

In a number of recent experiments the degree of spin 
polarization of quantum Hall fluid has been explicitly 
measured as a function of the ratio of Zeeman to Coulomb 



energy for a variety of fixed filling factors i^"— In this pa- 
per, we will focus on Refs.y and0, where detailed com- 
parison to our theoretical work is possible. 

In Ref. y a non-trivial net spin polarization is ob- 
served at the following filling factors in the lowest Landau 
level (LLL): u = 2/3, 3/5, 4/7, 2/5, 3/7 and 4/9. In these 
experiments the aim is to keep the filling factor fixed 
while varying the ratio of the Zeeman to the Coulomb 
energy. This is achieved by varying both the applied field 
strength and the density of electrons in the sample, keep- 
ing their ratio (hence v) fixed. For each filling factor, as 
the field strength is varied, the experiments report a set 
of plateau with constant spin polarization, punctuated 
by a series of relatively sharp transitions. Similar con- 
clusions can be inferred from data presented in Ref. |7| at 
filling factors i^ = 4/3, 7/5, 10/7, 8/5, 11/7 and 14/9. For 
high field and high electron density, the spin is polarized, 
but as the field strength and electron density are reduced, 
there are transitions to FQH states of successively smaller 
net spin polarization until, finally, some lower bound of 
spin polarization is reached, with the particular lower 
bound depending on the filling factor. These transitions 
are observed to occur at some critical values of the ap- 
plied field B"'^ and there is, therefore, a corresponding 
critical Zeeman energy per electron, i?^'' = giJ-sB"'* . 

From a theoretical perspective, an exceptional phe- 
nomenological understanding of the fractional quantum 
Hall effect (FQHE) has been acquired via the concept 
of the composite fermion (CF)i^ The key notion of CF 
theory is that the problem of strongly interacting elec- 
trons in a perpendicular magnetic field can be mapped 
onto the problem of non-interacting CFs in an effective 
magnetic flcld. The direction of the effective magnetic 
field can be cither parallel or antiparallcl to the physi- 
cal magnetic field. The CFs can be pictured as electrons 
bound to a certain number, p, of magnetic flux quanta. 
In the effective magnetic field, there exist effective Lan- 
dau levels — these are completely analogous to the LLs 



that occur for noninteracting electrons in the presence of 
a magnetic field. The FQHE of electrons can be inter- 
preted as the integer quantum Hall effect of CFs with p 
fluxes attached and occupying a certain number n of the 
effective LLs; this idea has proved to be very successful 
(see, for example, J. K. Jain's book, Ref. Q)- The prin- 
cipal set of filling factors encompassed by CF theory are 
given by 



V = n/(2pri ± 1) 



(1) 



where the sign here indicates the direction of the effec- 
tive field relative to the real magnetic field. Taking into 
account the two spin species, particle-hole conjugate ver- 
sions of these states occur at filling factors 2 — v. 

In CF theory, trial quantum Hall states with a spin 
degree of freedom can be constructed by simply associat- 
ing a spin degree of freedom with the CFs themselvesi^ 
CFs of each spin species can independently occupy a 
non-negative integer number n^ and n^ of effective LLs; 
the filling factor remains as in Eq. ([T|), but now with 
n ~ n^ + JT-j,. Consequently a whole scries of CF wave- 
functions is possible at each filling factor (a visualization 
of such a series of states can be found in Refs. [^ or[^ and 
a modified version is presented here in Fig. [T|). 

The experiments of Ref. y examined thep = 1 series for 
n = 2, 3, 4 with both positive and negative effective mag- 
netic fields: ly = 2/5, 3/7 and 4/9 in the positive effective 
magnetic field case and ly — 2/3, 3/5 and 4/7 for the nega- 
tive effective field case. Rcf. |7J examined the particle-hole 
conjugates of these states at filling factors 2 — i/. In each 
case, qualitatively at least, the predictions of CF theory 
appear to support experimental observations. We shall 
elaborate more on these qualitative predictions in Sec. 

El 

Transitions between CF ground-states with different 
spin polarizations can occur when the difference in 
Coulomb energy per electron between the two ground- 
state CF configurations compensates for the increase in 
Zeeman energy per electron due to the spin depolariza- 
tion. The differences in the Coulomb energies per elec- 
tron of various spin CF states at the same filling fac- 
tor thus can be related to the critical Zeeman energy 
per electron for transitions between the various spin- 
polarizations of the quantum Hall fluid. Calculating the 
Coulomb energy per electron for composite fermion trial 
wavefunctions. Park and Jain^ were then able to directly 
compare the experimental measurements and theoretical 
predictions of the critical Zeeman energies per electron 
for filling factors v = 2/5, 3/7 and 4/9 (the positive ef- 
fective field case). Their results for the predicted values 
of the critical Zeeman energy per electron for these three 
filling factors, for the most part, agree well with the val- 
ues measured in both Refs. |7| and y. These authors did 
not perform the same calculation and comparison for fill- 
ing factors v = 2/3, 3/5 and 4/7 (the negative effective 
field case). In this work we will present the results of 
our calculations for the negative effective field case, we 
shall evaluate the critical Zeeman energies for transitions 



between different spin states and then we will compare 
our results to the relevant experimental measurements. 

CFs have been the subject of a great deal of inves- 
tigation and there are now well-developed techniques 
which have been established for calculating the associ- 
ated Coulomb energy using the well-known Metropolis 
Monte Carlo procedure (sec. e.g., Refs. [9l- [lll) . Never- 
theless, for the series of states of interest here, with the 
exception of the spin polarized cases (see Ref. [ill ), the 
Coulomb energies have not yet been calculated. For suf- 
ficiently large numbers of particles the numerical evalu- 
ation of the trial wavefunctions with negative effective 
field turns out to be highly non-trivial — and much more 
complex than the positive effective field case. We have 
constructed an efficient new numerical algorithm to han- 
dle this situation. The details of the algorithm arc dis- 
cussed in Appendix iDl 

In this work, we calculate the interaction energy as- 
sociated with the CF trial wavefunctions using a simple 
Coulomb interaction potential, what we have called the 
Coulomb energy. Such an interaction would apply to 
a perfectly 2D geometry; however, a laboratory quan- 
tum Hall system cannot be considered perfectly 2D and 
more realistic model interactions must take into account 
the finite extent of the system in the direction perpen- 
dicular to the 2D plane. To study the effect of such a 
modification to the theory, we implement an interaction 
which takes into account finite-thickness effects (see, e.g., 
Ref. \12t \ . Finally, we calculate the Coulomb energy ap- 
propriate for the 2nd LL analogies of the LLL CF trial 
wavefunctions discussed previously (the 2nd LL being the 
next LL above the LLL i.e., filling factors v — v lll + 2). 
Using these results we make an additional prediction that 
for CF trial ground-states it is not energetically favorable 
to have non-polarized states in the 2nd LL with negative 
effective magnetic field. 

The structure of this paper is as follows: in Sec. Ullwe 
shall briefly summarize the qualitative predictions of CF 
theory and then we shall write down the explicit forms of 
the CF wavefunctions which we are interested in. In Sec. 
mil we shall present the results for the Coulomb energy of 
the various ground-state trial wavefunctions in the LLL 
and 2nd LL, the critical Zeeman energy predicted by CF 
theory for the LLL and 2nd LL and the results of the 
flnite thickness correction. We shall also compare our 
results to the experimentally measured values of the crit- 
ical Zeeman energy in the LLL. Finally, we shall make 
some remarks on our flndings in Sec. IIVI 



II. THEORY OF COMPOSITE FERMIONS 
WITH SPIN 

A. Qualitative predictions of CF theory 

Before we describe the explicit forms of the CF trial 
wavefunctions, we shall flrst briefly summarize the qual- 
itative predictions made by CF theory. The filling fac- 
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FIG. 1: A Summary of the qualitative predietions of CF theory (modified from a similar figure in Ref. |^. The figure 

shows the filling of effective LLs with spin-up and spin-down composite fermions. Different cases are labelled by 

their quantum numbers n-^ and nj^, the number of filled spin-up and spin-down LLs, e.g., (1:1) denotes n-^ = 1 and 

nj, = 1. The degree of spin polarization 7e is calculated using Eq. ^. 



tors of interest correspond to CFs with p = 1 fiux at- 
tached and occupying n = 2-4 effective LLs. In the 
case of negative effective field, using Eq. ((T|), these are 
v = 2/3,3/5,4/7, respectively, and in the case of posi- 
tive effective field these are ly = 2/5, 3/7, and 4/9, re- 
spectively. For each filling factor the set of possible 
spin-dependent states is deduced by considering all pos- 
sible non-negative integer values of n-^ and n^ satisfying 
n = n-^ + n^. For a system of N-f spin- up electrons and Ni 
spin-down electrons, we shall define the "degree of spin 
polarization" by 7e = rJ •rJ' ^-nd in the thermodynamic 
limit, where each effective LL contains the same number 
of electrons, it can be shown that (see Appendix [B]) 



7e = 



n^ + ni' 



(2) 



A description of the possible ground-states is summa- 
rized in Fig. [1] If the Zeeman energy is sufficiently large, 
then it is expected that the spin will be fully polarized, 
and so in high magnetic fields the system is pictured by 

I 



the rightmost diagrams in Fig. [T] We then reduce the 
applied field, keeping the filling factor fixed by lowering 
the electrons density at the same rate. When the crit- 
ical field i?"" is reached, then it is energetically favor- 
able for a transition to one of the states with lower net 
spin polarization, pictured in the diagrams successively 
to the left in Fig. [TJ This is energetically favorable as 
long as the difference in Coulomb energy between the two 
ground-state configurations compensates for the increase 
in Zeeman energy due to the spin depolarization. 

Comparing these qualitative predictions to the exper- 
imental results presented in Ref. y, we find that the pre- 
dictions for the number of transitions and for the degrees 
of spin polarizations are broadly correct as a first ap- 
proximation. In practice, the transitions are somewhat 
broadened and it is apparent that there are some small 
second-order plateau occurring between the main tran- 
sitions. These effects are not well understood and so 
presently we shall only focus on the leading-order effects. 
Further, we shall assume that the experiment is an ob- 
servation of ground-state quantum Hall behavior. 



B. Trial Wavefunctions for the LLL 



inates boundary effects. We have chosen to use the spher- 



In order to study the bulk properties of a quantum 
Hall ground-state we must choose a geometry which elim- 



ical geometry for this purpose, i.e., we shall study various 
finite-sized FQH states existing on the surface of a sphere. 
To effect a magnetic field perpendicular to the surface of 
a sphere we must place a magnetic monopole at its cen- 
ter such that the electrons see 7V$ flux quanta. As we 
increase the system size the total flux increases, but to 
fix the flux density at the surface of the sphere, the ra- 
dius of the sphere must also increase. On extrapolation 
to the thermodynamic limit, therefore, a flat geometry is 
recovered and the edge effects are eliminated. 

We shall now describe how CF states can be expressed 
as co-ordinate wavefunctions in the sphere geometry, in 
terms of the spinor co-ordinates Ui and Vi . (See Ap- 
pendix 12] for a more detailed discussion of the sphere 
geometry.) 



Spin polarized CFs 



The corresponding effective magnetic field thus can be 
aligned either parallel or antiparallel to the real mag- 
netic field. The intuition is that $„ represents an inte- 
ger quantum Hall wavefunction for composite fermions 
in an effective magnetic field and with effective filling 
n ~ ImiN^oo N/{2Q), i.e., the number of occupied ef- 
fective LLs is an integer n. For a finite sized system, 
N/(2Q) may be slightly shifted from its thermodynamic 
value. CF wavefunctions are often simply denoted as 
^PCF„ or ^PCF_„ where the sign corresponds to the sign 
of the effective field and n is now a positive integer^ We 
shall adopt this nomenclature for the remainder of our 
discussion. The CF states occur at filling factors given 
in Eq. ((!]). We shall always consider non-interacting CFs. 
Practically, the CF states are constructed as follows: 
for a system of N spin polarized electrons in the spher- 
ical geometry filling n CF LLs, the effective monopole 
strength is given by 



The CF wavefunction describing interacting electrons 
in magnetic field B can be succinctly expressed as 



^: 



p,7l 



Plll 



K^^n 



(3) 



where 



$0 = J]^ {UiVj 
i<j 



ijV,) 



and where -Plll denotes projection onto the lowest Lan- 
dau level (LLL). $„ is a Slater determinant of non- 
interacting single-fermion wavefunctions with an effec- 
tive magnetic fiux which we denote 2Q (i.e., the non- 
interacting fermions see 2Q flux quanta). The quantity 
Q is known as the effective monopole strength: it can 
take positive or negative integer or half-integer values.— 



±- 



iV- 



2n 



(4) 



with the sign depending on the sign of B^a. The mag- 
netic flux experienced by the electrons due to the mag- 
netic fleld B is then given by iV$ = 2p{N — 1) -I- 2Q [from 
Eq. ^]. In the sphere geometry, the single- particle CF 
eigenfunctions are the set of what are called monopole 
harmonics^ These monopole harmonics are eigenfunc- 
tions of the effective LL with eigenvalue n' = 0, ....n — 1, 
of the orbital angular momentum with eigenvalue I = 
\Q\, \Q\ + 1, ...IQl + n' , and of the z component of orbital 
angular momentum with eigenvalue m — —I, —l-\-l^ ...^l. 
It is simple to check that given Eq. (|3]) the total number 

of single particle states is J27=o (^(IQI +*) + !) = ^■ 
For Q < 0, the monopole harmonics are of the following 
form^^: 



^^'^.^ K,«.) = (-l)"'MQ^„,,„«)-Q+"«)-«-'"f](-l)-' (^ '^' ) (^ \Q\+nits ) «"0^««0"'-^ 

where i indicates the particle number and will run from I to the total number of particles A^. Here Mq^n' .m is an 
unimportant normalization factor. We can then write the Slater determinant as 



$. 



det 



Yl,n{y^^^^^) 



Yr 



0_|Q|("l.«l) 



rQ 



V-|Q|("JV,t'Jv) 



Y, 



I'^-IQI-i^"!'"!) 



Y^.qAun^vn) 



-Q 



^1 _|Q|_l(wAr,f'Af) 



^ 1 1^1 ( i(ui,Wl) ... y 1 1^1 , Jun.vn) 
I 



In order to implement the LLL projection required by Eq. ([3]), we follow the method introduced by Jain and 



KamillapiS, and then extended for the case of a negative 
effective field by Moller and Simonii. From a computa- 
tional perspective, for large system sizes, the LLL projec- 
tion of the Slater determinant is completely impractical 
(see Ref. [^ for details of the LLL projection). The is- 
sue can be circumvented by moving the Jastrow factor 
inside the determinant and then applying Plll to each 
of the resulting matrix elements first, before calculating 
the determinant: 



V'; 



p,n 



det 



PLLL{y^,mi^^,v,)J! 



(5) 



where 



n ^'^'^^ 



UjVt) . 



Although the result of this procedure is not mathemat- 
ically identical to Eq. ([3]), the resulting trial wavefunction 
would nevertheless describe fermions in the LLL, and the 
two prescriptions have been found to be extremely simi- 
lar in cases where they can be compared. An expression 
for 



Y, 



Q 



iu,,v,)J,^ PLLL{Yl%{u^,Vi)J!) 



in negative effective field was derived in Ref. [Ill and is 
given by 



t?™(^..«o«E(-i)' 



s=0 



2\Q\+n' 
\Q\+m + s 



d 
dui 



\Q\+m+s 



_d_ 

dvi 



\Q\—m+n ~s 



r 



Using this result, Eq. (O for the CF wavefunction can 
then be rewritten as 



2PCF_„ = ^p,„ = det \y^,,„,{u,, v,)Jf 



(6) 



In Appendix |D] we discuss an efficient technique for nu- 
merically evaluating such wavefunctions. 



CFs with spin 

If the CFs have a spin degree of freedom, then we must 
associate a spin degree of freedom with the single-particle 
monopole harmonic wavefunctions. Let us say we have 
N CFs with one of two possible spin polarizations, that 
is, N^ spin-up CFs and N^ spin-down CFs. These CFs 
can then independently occupy a number n-|- and n^ of 
spin-up or spin-down effective LLs. The value of p is 
independent of the spin degree of freedom, since it is 
not involved in the single particle monopole harmonic 
functions. The spin CF wavefunctions are of the general 
form: 



^p,(ri 



^LLL 



'I'^^*. 



,$. 



(7) 



where the Slater determinants $„^ and $„^ are formed 
from monopole harmonics with the following structures: 

The effective monopole strength is now given by: 



Q = ±- 



N^-72 



=:±- 



N_i-n 



'i 



2nt 2n_i. ■ ^ ' 

The possible eigenvalues of the monopole harmonics are 



0, 



■•,"-t' ^t 



101, -IQI 



and 771,-f- 



— /-|-, ■■■If, and similarly for the spin-down versions. Once 
we construct a Slater determinant of such states we can 
factor out the antisymmetric spin part of the wavefunc- 
tion, and we only need to specify the spatial parli^. 

In accordance with Jain's notation;^ we denote the se- 
ries of spin un-polarized CF wavefunctions by ^^CF(-„ „ ) 
or by ^^CF(_„ „ ) if the effective field is antiparallel to 
the magnetic field, with rif and ni being positive integers 
(it is not possible to have positive B^ff for one spin species 
and negative B^ff for the other) . The filling factor of the 
spin-dependent CF states is again given by Eq. ([Ijbut 
now with n = n-f -|- nj,. 

The final step is to project the wavefunctions onto 
the LLL; the final form of the spatial part of the spin- 
dependent CF wavefunction is, thus. 



^PCF(_„^,_„j = V'p.K.nj = dot [i;^,,„^(u,,w,)J«J det [^„'?^„^(m», w») J, 



(9) 



where the i index runs from 1 to N-f in the first determi- nant and from N-f + 1 to N^ in the secondJ^ Note that 



the Ji function is exactly as in Eq. © (i.e., it is a prod- 
uct over all spin-up and spin-down particles) and so in 
Eq. ([9]) each matrix element in each determinant depends 
on the co-ordinates of all of the particles. 



III. RESULTS OF MONTE CARLO 
SIMULATIONS 

In this section we shall present the results of our 
calculations of the Coulomb energy per electron asso- 
ciated with the seven principle spin CF trial ground- 
state wavefunctions in negative effective field: ^CF_4, 

2/^1? 2/^17 2/^17' 2/-^T7' 2/^17' 

^^(-3,-1)) '-'^(-2,-2)1 L^-f^_3, Oi^(_2,_i), Oi^_2, 

and ^CF(_i__i) . Using these results we shall deduce some 
quantitative predictions of CF theory for the critical Zee- 
man energy per electron at which transitions between 
different spin states occur. 



A. Coulomb Energy 

Our first calculation is the Coulomb energy per elec- 
tron, Ec, of the CF trial ground-state wavefunctions. In 
general, the Coulomb energy associated with a state de- 
scribed by a wavcfunction tp is calculated using 



(^|y|^) 



dvi...dYN\il)\ V{yi,...vn)-, 



(10) 



where V is the Coulomb potential (in units of e- / eIq): 



V 



N ^ 

^.R^. 



K, 



i<j 



Here Rij is the distance between pairs of particles in 
the sphere geometry. Throughout this paper we use the 
chord distance convention, Rf*'""^ = 2Rs\uiVj — VjUi\, 
where Rs is the radius of the sphere. Vbg is the poten- 
tial due to a uniformly changed positive background (we 
need to put this in so the overall system is electrically 
neutral) , and it is given by the self-energy of a uniformly 
charged sphere of charge -l-A^e plus the electrostatic en- 
ergy between the electrons and the uniformly charged 
sphere. Using the chord distance measure, Vbg = ^jj^ 
(seeRef.[il). 

One can numerically evaluate integrals of the form 
given in Eq. (jlOp using a Metropolis Monte Carlo pro- 
cedure (sec Appendix [Q for a brief description of how we 
implemented the Metropolis algorithm). The principal 
difficulty in using the Metropolis algorithm for our cal- 
culation is that the value of the wavcfunction must be 
re-calculated for many millions of different sample par- 
ticle configurations and, for CF trial wavefunctions, this 
process can be very computationally demanding. In or- 
der to make the procedure viable we require an efficient 
algorithm with which to evaluate the wavefunctions. A 



previously described algorithroii to evaluate CF wave- 
functions for the negative effective field case (see Ap- 
pendix [D]) is most computationally efficient for larger 
values of the effective LL number n and very inefficient 
for the smallest values of n, e.g n = 2 or n-j- or nj^ = 1 
or 2. For the purposes of studying the spin CF states 
we have seen that we often need to consider small values 
of n-|- or n^ and, consequently, until now, accurate cal- 
culation of energies for the spin CF states has not been 
computationally feasible. We were able to design an al- 
ternative algorithm which is most efficient for smallest n 
and less efficient for larger values of n. We shall discuss 
the details of our algorithm in Appendix |D] Using this 
new algorithm enabled us to calculate all of the results 
presented in this paper. 

The results of our calculation of the Coulomb energy 
per electron of the seven CF trial ground-state wavefunc- 
tions of interest are presented in Table |lll graphs of the 
extrapolations to the thermodynamic limit are presented 
in Fig. H 



B. 2nd Landau Level Coulomb Interaction 

An interesting consideration is to model equivalent 
wavefunctions in the 2nd LL, that is, the analogous states 
occurring at filling factor u — i^lll + 2. One approach 
would be to construct explicit wavefunctions at this new 
filling factor, however, the CF wavefunctions in anything 
other than the LLL are difficult to evaluate. A more ef- 
ficient alternative is possible, however: the problem of 
electrons interacting via the Coulomb interaction in a 
higher Landau level is mathematically equivalent to the 
problem of electrons in the LLL interacting via an ef- 
fective potential V'^^ir). An appropriate form for the 
effective potential is derived in Ref. [l5|, although, strictly 
speaking, this result was derived only for the disc geome- 
try. However, since we will take the thermodynamic limit 
at the end of the calculation, we can still use the same 
potential in the sphere geometry since the disc and the 
sphere geometries are the same in the thermodynamic 
limit. The effective potential is given by 



F-(.) = l + ^, 



(11) 



1=0 



with r in units of the magnetic length. 

The values of the coefficients q appearing here 
are deduced by equating the Haldane pseudopotential 
coefficientsi^ of the effective potential in the LLL, namely 



y<=«:== 



22'»+im! 7o 



rdrV''\ry^e-'-^l^, 



to the pseudopotential coefficients of the Coulomb inter- 
action in the second LL, namely 

Vl = J^°° qdqViq) (^1 (y ) ) ^^ (<?') e"^' , 



where 



2nd LL are given alongside the LLL values in Table |TI1 



V{q) 



rdrV{r)Jo{qr). 



In these expressions, L^ are Laguerre polynomials and Jq 
are Bessel functions. It is sufficient to work with only the 
first seven coefficients and we have checked that the addi- 
tion of more coefficients does not change the result signif- 
icantly. In order to determine the coefficients cq to cq in 
Eq. (fTTj) we must equate F,^" to T^^ for m = I, 3, 5, ... , 13 
(since the pseudopotentials vanish for m = 0,2,4,... 
etc.). Appropriate values of Ci are given in Table|lpi and 
we note that these differ very slightly from the equiva- 
lent values presented in Ref. [l5|. The effective Coulomb 
energies due to the modified interaction potential for the 

I 



Coefficient 


Value 


Co 


-50.36597363 


Cl 


87.38179510 


C2 


-56.08455086 


C3 


17.76579124 


a 


-2.971636200 


C5 


0.2513169758 


C6 


-0.008434843187 



TABLE I: Values for the coefficients in Eq. ([TTl)- We note 
that these differ very slightly from the equivalent values 
presented in Ref. [la- 
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FIG. 2: Extrapolations to the thermodynamic limit for the Coulomb energies per electron of various CF trial ground- 
state wavefunctions in the LLL. Only the linear extrapolations in 1/N are shown. Error bars are smaller than the 
data markers and so are not plotted. 



System Size 


Af* 


Energy (LLL) 


Energy (2nd LL) 


20 


37 


-0.48453(2) 


0.5773(1) 


24 


44 


-0.48508(2) 


0.5844(1) 


28 


51 


-0.48548(2) 


0.58980(9) 


32 


58 


-0.48575(2) 


0.59386(9) 


36 


65 


-0.48614(4) 


0.5976(2) 


ooi 




-0.4880(1) 


0.6222(7) 


oo^ 




-0.4887(8) 


0.629(2) 


oo^ 




-0.4881(2) 


0.623(1) 


oo^ 




-0.4883(5) 


0.625(2) 



(a) 2Ci^_4 State 



System Size 


iV$ 


Energy (LLL) 


Energy (2nd LL) 


14 


25 


-0.49595(4) 


0.6166(2) 


18 


32 


-0.49483(3) 


0.6201(2) 


22 


39 


-0.49398(3) 


0.6224(1) 


26 


46 


-0.49353(3) 


0.6240(1) 


30 


53 


-0.49309(4) 


0.6246(2) 


34 


60 


-0.49277(4) 


0.6253(2) 


38 


67 


-0.49266(5) 


0.6263(2) 


oo^ 




-0.49062(6) 


0.6318(3) 


oo^ 




-0.4906(2) 


0.6305(9) 


oo--' 




-0.49059(9) 


0.6315(4) 


oo* 




-0.4907(1) 


0.6311(5) 



(b) ^CF, 



(-3,-1) 



State 



System Size 


iV$ 


Energy (LLL) 


Energy (2nd LL) 


12 


21 


-0.50276(3) 


0.6430(1) 


16 


28 


-0.49986(3) 


0.6401(1) 


20 


35 


-0.49809(3) 


0.6388(1) 


24 


42 


-0.49677(3) 


0.6375(1) 


28 


49 


-0.49597(5) 


0.6366(2) 


32 


56 


-0.49538(5) 


0.6365(2) 


36 


63 


-0.49488(4) 


0.6355(2) 


40 


70 


-0.49442(4) 


0.6354(2) 


oo^ 




-0.49088(5) 


0.6321(2) 


2 
CO 




-0.4908(2) 


0.6321(5) 


OO^ 




-0.49082(6) 


0.6321(2) 


oo* 




-0.4908(1) 


0.6319(3) 



(c) ^CF, 2,-2) State 



System Size 


Ni, 


Energy (LLL) 


Energy (2nd LL) 


12 


21 


-0.49843(3) 


0.6247(1) 


15 


26 


-0.49810(3) 


0.6325(1) 


18 


31 


-0.49801(3) 


0.6382(1) 


21 


36 


-0.49781(2) 


0.6420(1) 


24 


41 


-0.49765(2) 


0.6448(1) 


27 


46 


-0.49756(2) 


0.6472(1) 


30 


51 


-0.49754(4) 


0.6494(2) 


33 


56 


-0.49745(4) 


0.6509(2) 


36 


61 


-0.49741(4) 


0.6522(2) 


oo^ 




-0.49690(4) 


0.6657(2) 


oo^ 




-0.4968(1) 


0.6672(5) 


oo^ 




-0.49688(6) 


0.6660(2) 


oo^ 




-0.49679(5) 


0.6661(3) 



System Size 


iV$ 


Energy (LLL) 


Energy (2nd LL) 


8 


13 


-0.52139(5) 


0.7091(2) 


14 


23 


-0.51188(4) 


0.6930(2) 


20 


33 


-0.50821(4) 


0.6868(2) 


26 


43 


-0.50611(4) 


0.6833(2) 


32 


53 


-0.50488(4) 


0.6814(1) 


38 


63 


-0.50416(6) 


0.6805(2) 


44 


73 


-0.50328(6) 


0.6794(2) 


oo^ 




-0.49938(6) 


0.6724(2) 


oo' 




-0.4994(2) 


0.6733(4) 


oo^ 




-0.4994(1) 


0.6728(3) 


oo* 




-0.4994(2) 


0.6732(5) 



(e) 'CF,2, 



State 



(d) ^CF_3 State 



TABLE IL Coulomb energy per electron, Ec, for trial wavefunctions at filling factors v = 2/3, 3/5, and 4/7, for various 
system sizes; some of these data are plotted in Fig. [2l Results are given for the Coulomb interaction in the LLL, 
and for the modified Coulomb interaction in the 2nd LL. Energies are calculated using chord distance measure in the 
sphere geometry, as defined in the text, and are stated in units of e^/e/g. The extrapolations to the thermodynamic 
limit were calculated using three different methods, indicated by the superscript 1-4: 1 is with linear extrapolation in 
1/A^, weighted by the statistical errors on each data point; 2 is with quadratic extrapolation in 1/iV, weighted by the 
statistical errors on each data point; 3 is with a weighted linear extrapolation in 1/A^, excluding the smallest system 
size; and 4 is with a weighted linear extrapolation in 1/N excluding the two smallest system sizes. The different cases 
are calculated in order to estimate the possible error in the extrapolations. We note that for the ^CF-4, ^CF-3 and 
^CF_2 states, our numbers are in good agreement with those obtained in Ref. |Tl|. 



System Size 


iV$ 


Energy (LLL) 


Energy (2nd LL) 


6 


9 


-0.53899(4) 


0.7716(2) 


10 


15 


-0.53047(3) 


0.7677(2) 


14 


21 


-0.52675(3) 


0.7669(2) 


18 


27 


-0.52468(3) 


0.7667(1) 


22 


33 


-0.52342(2) 


0.7663(1) 


26 


39 


-0.52257(5) 


0.7665(2) 


30 


45 


-0.52181(5) 


0.7663(2) 


34 


51 


-0.52129(4) 


0.7663(2) 


38 


57 


-0.52096(4) 


0.7665(2) 


oo^ 




-0.51755(2) 


0.7649(3) 


oo^ 




-0.51751(5) 


0.7666(1) 


oo^ 




-0.51753(3) 


0.7657(1) 


oo^ 




-0.51754(5) 


0.7660(1) 



System Size 


iV* 


Energy (LLL) 


Energy (2nd LL) 


8 


11 


-0.55748(4) 


0.8754(2) 


12 


17 


-0.54506(4) 


0.8355(1) 


16 


23 


-0.53906(3) 


0.8172(1) 


20 


29 


-0.53536(7) 


0.8061(1) 


24 


35 


-0.53309(3) 


0.7991(1) 


28 


41 


-0.53133(5) 


0.7942(1) 


32 


47 


-0.53016(3) 


0.7904(1) 


36 


53 


-0.52914(3) 


0.7874(1) 


40 


59 


-0.52837(3) 


0.7850(1) 


oo^ 




-0.52102(8) 


0.7621(5) 


oo^ 




-0.52144(6) 


0.7652(3) 


oo^ 




-0.52118(5) 


0.7633(2) 


oo^ 




-0.52123(6) 


0.7636(2) 



(f) 2Ci^_2 State 



(g) ^CF( 



(-1,-1) 



State 



TABLE IL (Continued) 
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C. Comparison with Experiment 

At this point we are now ready to calculate the pre- 
dicted values for the critical Zeeman energy per electron 
E"^'*" for transitions and to make comparisons with the 
experiments. To study each transition, e.g., at filling 
i/=y,7e = l — > 2 (for which the CF theory prediction 
comes from the transition ^CF-4 — s- ^Ci^(_3__i)), we 
need to compare the two key energy scales: first, there 
is the difference in Coulomb energy per electron, A.Ec, 
between the two participant ground-state trial wavefunc- 
tions; second, there is the difference in Zeeman energy 
per electron, AEz, which is due to the difference in the 
net spin polarization of the two participant ground-state 
trial wavefunctions. The condition for a transition is that 
AEc = AEz at the point of transition. 

The difference in Coulomb energy per electron, AEc, 
is calculated directly from our results for the Coulomb 
energy per electron associated with each trial ground- 
state wavefunction — these differences are listed in Table 
mil The difference in Zeeman energy per electron as we 
go through a transition is equal to the energy E"^^* to flip 
a single spin in a magnetic field i?'"" multiplied by the 
proportion of spins r which need to be flipped as we go 
through that transition. In the thermodynamic limit the 
proportion of flipped spins is r = 1/?t-, where n is the total 
number of filled effective LLs of the CF wavefunctions 
partaking in the transition (see Appendix [B]) . Thus, at 
the transition we have AEc = AEz = E"^'* /n and this 
leads to the relation 



Ez" = nAEc. 



(12) 



Using Eq. (jl2[) we have calculated predicted values for the 
critical Zeeman energy, and these predictions are listed 



in Table Uni 

Values of i?f " can also be deduced from the exper- 
imental results presented in Refs. and 0. In Ref. [j 
the values are derived from the measurements of the de- 
gree of spin polarization je as a function of magnetic 
field strength at fixed filling factor. Accompanying each 
transition is a broadened step in the degree of spin po- 
larization. We take the critical field B"" to occur at the 
center of the step, and we take the experimental error 
to be roughly the half-width of the broadening (which 
corresponds approximately to an error of ± 0.001 in the 
Zeeman energy). In Ref. |7| we infer the critical fields 
for spin transitions at filling factors v = 2/3,3/5 and 
4/7 from the critical fields for spin transitions observed 
in particle-hole conjugate states occurring at filling fac- 
tors v = 4/3, 7/5 and 10/7. Particle-hole conjugation 
does not, in principle, affect the Coulomb or Zeeman en- 
ergy associated with the trial CF wavefunctions and so 
the predictions for the critical Zeeman energy are iden- 
tical. The aim of the experiments described in Ref. is 
to vary the ratio of Zeeman to Coulomb energy keeping 
the filling factor fixed. In order to achieve this the elec- 
tron density is kept constant and the applied field 5'°' 
is tilted by an angle 9 from the verticali^ and increased 
in magnitude simultaneously, thus fixing the component 
of the field perpendicular to the plane, B± = _B'°' cos6'. 
The signature for transitions between states of different 
spin polarization is a peak in the ratio of the longitudinal 
resistivity at a given filling factor to the longitudinal re- 
sistively at filling i/ — 3/2. We take the critical field B"" 
to occur at the center of the peaks and we take the ex- 
perimental error to be the half- width of the peaks (which 
corresponds approximately to an error of ±0.001 in the 
Zeeman energy) . Using these values for B""* we calculate 
the experimentally measured value of the critical Zeeman 
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energy via E"^"' = gfisB""^, and then convert to units of 
e^/elo for comparison with our theoretically derived val- 
ues in Tabic IIIII (for this calculation we take the g factor 
of GaAs to be 5 = —0.44 and the relative permittivity to 

I 



be e = 12.6). 

In Fig. [3] we have plotted the predicted and measured 
values of i<Jf " given in Tabic IIIII as a function of the 
parameter n. 



Filling 


Transition 


CF Theory Prediction 


AEc 


Et' (Predicted) 


Ff " (Kukushkin et al.) 


E°£* (Du et al.) 


4 
7 
4 
7 
3 
5 
2 
3 


-v - 1 ^ 1 

7. = i ^ 

7e = 1 ^ 


CF(_3__i) — >■ CF{_2,-2) 

^CF-2 -^ 'CF(_i._i) 


0.0027(2) 
0.0003(1) 
0.0025(1) 
0.00363(7) 


0.0107(8) 
0.0010(4) 
0.0074(3) 
0.0073(1) 


0.013(1) 
0.007(1) 
0.012(1) 
0.0088(1) 


0.018(1) 
0.010(1) 
0.017(1) 
0.015(1) 



TABLE III: For each filling factor in the negative effective field case, the table shows the possible spin-transitions 
labelled by their degree of spin polarization 7e [defined in Eq. ^], the difference in Coulomb energy per electron 
between the two possible ground-state configurations, and the corresponding prediction for the critical Zeeman energy 
per electron calculated using Eq. (|12[). We use the extrapolation scheme that gives minimum uncertainty in the 
extrapolated value. The relevant experimentally derived values for the critical Zeeman energy taken from Kukushkin 
et al. (Ref. y) and from Du et al. (Ref. 0). We explain in the text how the experimental values are deduced. All 
values are given in units of e"^ /eIq. 
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FIG. 3: Predicted and measured values for the critical Zeeman energy per electron from Table lllll plotted against the 
reciprocal number of filled effective LLs, 1/n. Transitions arc labelled by their degree of spin polarization 7e [defined 
in Eq. ([2])] . In the text we explain how the theoretical predictions of the Critical Zeeman per electron are calculated. 
Experimental values are taken from Kukushkin et al. (Ref. y) and from Du et al. (Ref. 0) . We explain in the text 
how the experimental values are deduced. Circled sets of points correspond to the same transition. For the series 



of ^CF_ 



'CF, 



_i) transitions (i.e., those transitions for which higher value of 7e is 1) there is a trend for 



the critical Zeeman energy to increase as n increases, 
energy. 



In all cases theory tends to underestimate the critical Zeeman 
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D. Finite Thickness Correction 

In the experimental investigations of the FQHE dis- 
cussed in Refs. y and 0, the 2D geometry is reahzed in 
a GaAs-Al3;Gai_2,As heterojunction setup. Due to the 
finite width of the potential well in hetcrojunctions, the 
geometry cannot be considered perfectly 2D and real- 
istic wavefunctions must have some finite extent in the 
direction perpendicular to the 2D plane. The impact 
of the finite thickness correction can be calculated sim- 
ply by modifying the interaction potential with which 
we evaluate the ground-state energy. Appropriate mod- 
ified potentials are discussed, for example in Ref. [Tj. 
Here we implement the following potential, known as 
the Fang- Howard potential, which is applicable to GaAs- 
Ala:Gai_2;As heterojunctions. Note that this result was 
derived only for the disc geometry, however in the ther- 
modynamic limit, the result will be correct for the sphere 
geometry also. 



KH(r,d)=V;V.+ / dkV^^{k,d)Mkr)k, (13) 

"'0 



with 



V^^ik,d) 



e\ 9 24: + 9kd + {kdY' 
~8fc (3 + kdf 



The potential is a function of a thickness parameter d/lo, 
which characterizes the extent of the wavefunction in the 
dimension perpendicular to the plane; d = corresponds 
to a perfect 2D geometry. Note that the value of the 
background contribution to the potential, V^^, will de- 
pend on this new potential and will not be the same as 
for the Coulomb potential case. Using this modified in- 
teraction potential we have calculated the ground-state 
energy in the thermodynamic limit for a number of dif- 
ferent values of d. Graphs showing the thermodynamic 
extrapolation of the ground-state energy as a function of 
the thickness parameter are presented in Fig. 21 

The value of d for a particular GaAs-Ala;Gai_^As het- 
erojunction can be estimated using known data about the 
system (see Ref. [20) : 



1 /487rm^eV* 

3 V Kscft^ 



(14) 



where uiz is the effective mass of the electron in the GaAs 
in the direction perpendicular to the plane of the hetero- 



junction, which is 0.063me; ksc is the static dielectric 
constant of GaAs, which is 12.6 and p* for an un-doped 
system is given by 32 Ps, where ps is the areal density of 
electrons in the inversion layer. The value must then be 
expressed in units of the magnetic length ^o- 

In the experiments of Kukushkin et al. (Ref. |3|) the 
system is setup such that the ratio of the ps to B"'* is 
fixed by the filling factor, i.e., ps — vB""^ /4)q (here 0o 
is the flux quantum ^o = hc/e), and, thus, in Eq. P^ 
the value of d is a function of B"'^ . An additional depen- 
dence of d on B°"^ enters through the magnetic length. 
Theoretical predictions for B"" can be derived from the 
values of the critical Zeeman energy given in Table IIIII 
or, if we consider a system with non-zero d, from the 
results presented in Fig. |4l Clearly, then, the theoreti- 
cal prediction for _B"" itself depends on the chosen value 
of d, estimated from Eq. P^ . For each type of tran- 
sition we can determine theoretical predictions for the 
values of d by requiring that the above conditions are 
self-consistently satisfied. We have tabulated these self- 
consistent estimates in Table ITVl 

In the experiments of Du et al. (Ref. 0) there are 
two slight modifications to our estimation of d. First, 
in these tilted field experiments the value of ps is fixed 
at 1.13 X 10^^ cm~^. Second, the effective magnetic 
length due to the in-plane field _B'°' places an additional 
limit on the extent of the wavefunction perpendicular to 
the plane. The in-plane field can be determined by the 
requirement that the filling factor is fixed and the num- 
ber density of electrons is known, i.e., we require that 
i3_L = Ps^ajv; the value of the total field is derived from 
the prediction of the critical Z eeman energy; th e in-plane 
field B\\ is given by B\\ = ^ {B'°'f -{Bx_f. In fact, 
since our predictions for the critical Zeeman energy are 
relatively small compared to the values observed in Du's 
experiment, we find that our predicted values for 5'="* 
are smaller than the values B^ consistent with the ex- 
perimental electron density and filling factor. It is, there- 
fore, not possible in this case to construct a self-consistent 
value for Bn . We can, however, use the experimentally 
observed values of the critical field to make an estimate 
of the value of d in units of Iq for each transition using 
Eq. HID. 

In Table IIVI we present our estimates of d appropriate 
for each experimental setup described in the preceding 
two paragraphs and we present the corresponding theo- 
retical predictions for the critical Zeeman energy per elec- 
tron calculated using the potential described in Eq. (|13D . 
In Fig. [5]we plot the modified results for the critical Zee- 
man energy per electron given in Table IIVI as a function 
of n as in Fig. [3] 
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Thickness Parameter d [1 ] 

(a) Interaction energy per electron associated with a selection of the CF ground-states in the 
thermodynamic limit due to a modified Coulomb potential Vfh ('", t^) given in Eq. (|13p , plotted 
as a continuous function of the thickness parameter d. The associated error bars for these 
extrapolations are not drawn as they are too small to see on the plot. 



0.012 




Thickness Parameter d [1 ] 

(b) Theoretical predictions for the critical Zeeman energy per electron due to a finite-thickness 
potential Vfh ('', d) given in Eq. (J13p , plotted as a continuous function of the thickness parameter 
d for selected CF states. The critical Zeeman energy is calculated using Eq. (fT^ . Representative 
error bars for a selected subset of the data are drawn at intervals on the plot. Note that the 
data have been smoothed between d — 0.2 and 0.3 in order to remove a numerical artifact. 



FIG. 4: Graphs showing the impact of the finite thickness correction on our results. The energy is given in units 
of e^/e/o and the thickness parameter is given in units of Iq. To obtain these curves, weighted extrapolations to the 
thermodynamic limit are calculated for the sets of data points at finite N for 501 discrete values of the parameter 
d between and 5 and then the curves are interpolated. We use the extrapolation scheme that gives minimum 
uncertainty in the extrapolated value (i.e., possibly with the smallest system size removed from the extrapolation). 
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Kukushkin ei al. 


Du et al. 


Filling 


Transition 


CF Theory Prediction 


d 


£,crit 


d 


ET' 


4 
4 

3 

5 
2 
3 


7e = 1 ^ i 
7e = ^ ^> 
7e = 1 ^ i 
7e = 1 ^ 


^CF.4 ^ 'C7F(_3,_i) 

C-P(-3,-l) — >■ C'F(_2,_2) 

^CF.3 ^ 'C7F(_2,_i) 
^CF-2 -^ 'CF(_i,_i) 


1.42 
0.865 
1.33 
1.44 


0.0070(4) 
0.0016(3) 
0.0060(2) 
0.0063(1) 


0.595 
1.20 
0.636 
0.776 


0.0096(6) 
0.0014(3) 
0.0076(2) 
0.0078(1) 



TABLE IV: The table presents estimates of the finite thickness parameter d in units of /o for different transitions and 
filling factors in negative effective field. Transitions are labelled by their degree of spin polarization 7^ [defined in 
Eq. ([2])] . The estimates arc calculated using Eq. ([Ti|) and experimental data taken from Kukushkin et al. (Ref. y) or 
from Du et al. (Ref. [3). The table also presents modified theoretical predictions for the critical Zecman energy per 
electron in units of e^/e/o calculated by using the estimates for d with the modified potential given in Eq. [T31 
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FIG. 5: Predicted values for the critical Zeeman per electron modified by the inclusion of finite thickness effects 
(data from Table lIVp and measured values for the critical Zeeman energy per electron from Table |IlT] plotted against 
the reciprocal number of filled effective LLs, 1/n. Transitions are labelled by their degree of spin polarization 7e 
[defined in Eq. ([2])]. In the text we explain how the theoretical predictions of the critical Zeeman energy per electron 
are calculated, taking into account a finite thickness correction to the potential appropriate to the conditions of the 
different experimental systems considered. Experimental values are taken from Kukushkin et al. (Ref. |3|) and from 
Du et al. (Ref. 0). We explain in the text how the experimental values are deduced. Circled sets of points correspond 
to the same transition. As in Fig. |31 the theory somewhat underestimates the critical Zeeman energy. 
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IV. CONCLUSIONS 



Looking at Table lllli we conclude that the predic- 
tions made by CF theory for the critical Zeeman en- 
ergy agree moderately well with the experimental values: 



we find a good agreement for the ^CF_2 to ^Ci^(_i,_i) 
transition, but we find that the predictions for the 



Ci^_3 to ^CF,2 -1), the ^CF_i to ^CF, 



(-3,-1) 



and the 



Ci^(_3__i) to Ci^(_2,-2) transitions are, respectively, at 
least a factor of 1.6, 1.2 and 7 away from the experimen- 
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tal values. Looking at Fig. |31 wc notice a trend in both 
sets of experimental data for the critical Zeeman energy 
to increase with increasing n (excluding the ^CF(_3__i) 
to ^Ci^(_2,_2) transition). This trend is also present in 
the theoretical predictions of CF theory. Compared with 
the analogous set of results for the positive effective field 
case (fining factors v = 2/^,'i/7 and 4/9, see Ref. Q), 
those results tend to agree much better with experiment 
for the series of ^CFn — > '^CF(^n-i,iyi however, there re- 
mains a similar large discrepancy with experiment for the 
^CF(3^i) to ^CF(2,2) transition at filling 4/9. 

The potential impact of the finite thickness correction, 
for small values of the thickness parameter {d < 0.5), is to 
increase the differences in the predicted critical Zeeman 
energy per particle between the adjacent states typically 
by 20-30%. For larger values of the thickness parame- 
ter (d > 0.5 ) the predicted critical Zeeman energy per 
particle decreases. Based on our estimates of the value 
of d, given in Table IIVI we conclude the the effect of 
the finite thickness correction would in fact be to lower 
our predicted values for the critical Zeeman energy per 
electron. We do however, correctly predict that the crit- 
ical Zeeman energy should be greater under the experi- 
mental conditions of Du et aVs setup compared to the 
conditions of Kukushkin et aVs setup, with the excep- 
tion of the ^CF(_4)-to-^CF(_3._i) transition (see Fig. [5]). 
Nevertheless, the finite thickness correction clearly does 
not account for the discrepancy between the theoretical 
predictions and the experimental measurements. Other 
factors which have been neglected are LL mixing, finite 
temperature effects, and the effect of sample impurities 
and disorder. Another source of error is that the pre- 
cise value for the Coulomb energy in the thermodynamic 
limit depends on how one does the extrapolation: linear 
and quadratic regression, or regression taking into ac- 
count the relative errors on each data point, or not. This 
issue is particularly acute for the n = 4 case because, ac- 
cording to Eq. ((T2|) . we must multiply the differences in 
Coulomb energy by a factor of 4 to obtain the prediction 
for the critical Zeeman energy. 

Turning to the results for the analogous states occur- 
ring in the 2nd LL, we can deduce an interesting pre- 
diction: it would appear that for n = 4 and n ~ 3 the 
fully polarized states i^CF^^ and ^CF^^, respectively) 
lie lower in energy that any of the non spin-polarized 
states with the same filling factor, although for the n = 2 
case the results are inconclusive. For the states at filling 
factor 2 -I- 3/5 and 2 + 4/7 and perhaps for the states at 
filling factor 2 -I- 2/3 as well, the prediction of CF theory, 
therefore, would be that it would not be energetically 
favorable to depolarize the spin as the field strength is 
reduced from its highest value and we would expect to 
observe no spin transitions at all in these cases. An al- 
ternative possibility, which is particularly likely for the 
2 -|- 2/3 case, is that the energy differences are very small, 
and that would imply a very low value of the critical field 
as compared with the analogous cases in the LLL. An 
experimental observation of extensive non-polarized be- 



havior at filling factors 2 -I- 2/3, 2 -I- 3/5 or 2 -I- 4/7 would 
suggest the existence of ground-state wavefunctions not 
predicted by CF theory. 

The most recent experimental studies of the spin po- 
larization of FQHE states have focused on the question 
of the spin polarization of the u = 5/2 statei^i^i From 
the perspective of CF theory, the 5/2 sate is described by 
CFs with zero effective magnetic field^. Given our above 
conclusion that 2nd LL states described by CF theory 
would tend to favor a spin polarized configuration, we 
speculate that a spin polarized configuration may also 
be favorable at 5/2. This would be in agreement with 
exact diagonalizations of smaller systems^^i^i. 

Note added. As this manuscript was being prepared 
for publication we became aware of some very recent ex- 
periments to determine the spin polarization of quantum 
Hall states in the 2nd LL at filling factor 2 + 2/3, (see 
Ref. l28r ). In that work, evidence is presented indicating 
that a spin transition of the form we have been discussing 
here does occur at filling 2 + 2/3. The result seems to 
agree with the predictions of CF theory that we have 
enumerated here; however, it should also be pointed out 
that the ground-state wavefunctions may not necessarily 
be of the CF type (see Ref. HI" 
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Appendix A: Geometry 

The FQHE occurs in 2D systems in a perpendicular 
magnetic field. Real systems are, of course, finite in size 
and, thus will have edges; however, for our investigation, 
we are concerned only with the bulk properties of the 
ground-state and do not want to take into account edge 
effects. One method to eliminate edge effects, from a 
theoretical perspective^^ is to place the quantum Hall 
system on the surface of a hypothetical sphere of radius 
Rs- The surface of a sphere is described by a spherical 
co-ordinate system fl = {9,(j)) with a fixed radius Rs- 
Out of convenience we choose to write wavefunctions in 
this geometry using a pair of complex spinor co-ordinates 
u, V such that 



ict>/2 



i4,/2 



A magnetic field B perpendicular to the surface of the 
sphere can be realized by placing a magnetic monopole 
at the center of the sphere: the total magnetic flux is 
N^(j)Q = AttR^B (here (f>Q is the flux quantum 0o = he/ e). 
The radius of the sphere is then Rs = y^N^/2, in units 
of the magnetic length 



'-'5) 



1/2 
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The flux is related to the fiUing factor by A^$ = N/v — 
S, where N is the number of electrons and S is called 
the shift. For spin-polarized CF states we have already 
derived (in Sec. H]) the result 7V$ = 2p{N-l)+2Q, where 



Q 



± 



N- 



For the non-spin polarized CF states we 
, along with the definition N ^ N^ + N^ , 



2n 

can use Eq. 

to find an expression for Q in terms of N 



Q = ± 



2(nt + %) 



(Al) 



The total flux is then again given by iV$ = 2p{N—l)+2Q. 
In each case, in the thermodynamic limit, the sphere 
radius tends to infinity and so we recover a 2D plane ge- 
ometry but without edge effects. The length scale set by 
the magnetic length Iq is generally much smaller than the 
sphere radius, and so we argue that the physics remains 
independent of the system size and, therefore, taking the 
thermodynamic limit is valid. 



Appendix B: Degree of Spin Polarization and 
Proportion of Flipped Spins 

In this appendix we shall derive two results used in the 
paper: the degree of spin polarization 7 and the propor- 
tion of spins T which must flip when a spin transition 
takes place. 



1. Degree of Spin Polarization 

The degree of spin polarization is defined by 



7 



N^ + N^- 



Using Eq. ^ and the fact that N = N-^ + Ni wc can 
deduce for example that 



N-f-Ni = N - 



2n| 



714, + n^ 



{N + n^{ni -n^)). 



and so we have 



'H + "t \ ^ 

In the thermodynamic limit {N — > 00) we recover Eq. [51 

2. Proportion of Flipped Spins 

The proportion of electrons, r, occupying the effective 
LL labelled by index n' and with 2Q effective flux, is 
given by 

2{\Q\+n') + l 
^ = N • 



Substituting the expression for Q given in Eq. (jAip we 
have 




In the thermodynamic limit {N — ^ 00) this expression 
reduces to 

1 1 



nt + n; n 
which is independent of n' . 



Appendix C: Monte Carlo Algorithm for a Quantum 
Hall Fluid 



In this work we apply the Metropolis Monte Carlo 
algorithm^^ to evaluate the ground-state energy of var- 
ious trial wavefunctions. We wish to evaluate expecta- 
tion values of operators A with respect to the co-ordinate 
wavefunctions ^p (ri, ...tat): 



A)^ 



J dri...drNi/;* (ri, ...fat) AiIj (ri, ...tn) 

Jdri...drN |V^ (ri, ...rjv)|^ 

The Metropolis Monte Carlo procedure works by sta- 
tistically sampling configurations of the co-ordinates 
{ri, ...Tat} drawn from the probability distribution 



PAT (ri,...rAr) 



\ip{ri,...rN)\' 



Jdri...drN IV" (ri, •■•I'Ar)! 



2 • 



The expectation value of the operator is then estimated 
using A's sets of co-ordinate samples: 



1 



N, 



A) = -^ '^ip {ri, ...tn)* Alp {ri, ...tn) 



It can be shown that the standard deviation behaves as 
a/ < A >^ l/v^/Vg. In our simulations we typically used 



Appendix D: Algorithm for Numerical Evaluation of 

Composite Fermion States in Negative Effective 

Field 



Composite fermion wavefunctions have been intensely 
scrutinized using MC methods<^"— The principle diffi- 
culty for CF wavefunctions is the procedure for LLL pro- 
jection. The method for doing this projection for CF 
states in negative effective field was introduced in Ref. [ll|, 
and the resulting form of the generalized spherical har- 
monics is repeated here in Eq. ^. Once written in this 
form, the key difficulty lies in the evaluation of the mul- 
tiple derivatives of the Jastrow factor. 



_d_ 
duj 



\Q\+m+s 



_d_ 

dvi 



\Q\—in+n—s 



Jf 
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with 

Ji = J^ {UiVj - UjVi) . 

A procedure for evaluating these derivatives is given in 
Ref. [29 (see also Ref. Ig), and this method is also used in 
Ref. [ll|. Briefly, the method is as follows: first, we pull 
the Jastrow factor through the derivatives and write 



Jf 



f/' 



|+m+s-,~r|Q|— ni+n — s-i 



where 



f] ~ J-P— JP V- - rP— jP- 



we then introduce 

N 
k=l 



Vk 



UiVk - ViUk 



-Uk 



U,Vk ~ V,Uk 



from which one can deduce the recursion relations 

^/,(a, /3) = -(a + /?)/,(« + l,/3). 



Using these results, one can calculate a series of relations, 
for example. 



c/a = p/,(i,o) ,[/fi=pV«(i,o)2-p/,(2,o) ,. 



For CF wavefunctions in negative effective field we 
must take up to 2\Q\ + n ^ (■^) — 2n derivatives with re- 
spect to both Uj and Vj. As we take more derivatives, the 
results of this method become increasingly complicated, 
particularly for smaller values of n. 

We have determined an alternative method to evaluate 
the derivatives, which is at its most effective in precisely 
the regime where the current method runs into difficul- 
ties. We shall present our result for the case of p = 1 
only; however, higher p cases could be constructed along 
the same lines. 

Let us re-state exactly what we need to evaluate, leav- 
ing out the un-important constant factors as they can be 
absorbed into the normalization — for each element of a 
N-by-N Slater matrix we need to evaluate 



n 
s 



s=0 

It is insightful to expand out the J,; product into a sum as follows 

N N / 



2\Q\+n' \ , „'_, /^\ l'2l+"+-^ (_d_ ^ 101-™+"'-^ 
IQI+m + sP*''' \duj \dv, 



n ("*"j ~ "j"«) • 



j¥'i 



^ N \ /JV-1 \ 

u¥' J \t=o J 



where ej denotes the degree t elementary symmetric polynomial in the A^ — 1 variables Vj /uj for j ^ i, 

e\ = et,N-i{vi/ui,...Vj/uj...,VN/uN), 
for j ^ i. The elementary symmetric polynomials are defined by 



e„,jv(a;i,...,xjv) = { o<ji<j2<...<»™<w 





Xi-^...Xi^ m < N 



otherwise. 
With this form for Ji we now evaluate all of the necessary derivatives, which leaves us with: 

t?„(».,.w«Ei.(-i)-(";)(|^i«i::;)x 

N-l-{\Q\~m+n'-s) i, , w (JV-l-t)! N ^l-t-{\Q\-m) t\ 



l^t=\Q\+m+ 



% ( -\\t \^^ — 1 — r;i 1 

^t\^) N-\-t~{\Q\-m+n'~s)\'^i 



{t-{\Q\+m.+s)Y t 



t-(|Q|+m) 



(Dl) 



(D2) 



The expression in Eq. (|D2p appears complicated at first sight, but notice that the sum over t contains at most 
A' — 1 — (^) -I- 2n terms. If the CF state only fills a small number of effective LLs (i.e., n = 1, 2, ...), then the number 
of terms in the sum is actually quite small. 
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A key result is that the elementary symmetric polynomials can be calculated recursively using one of Newton's 
identities rl 



e-m,N 



_. m 

{Xi, ...,Xn) ^ —}{-iy'pr,N {xi,...,XN)em~r,N (xi, ...,Xn) : 



^N 



where Pr^N (a^i, ••■, xn) = X]i=i -^^ ^^'^ ^^'^ power-sum polynomials. Also, note the following recursive identity: 

em,N-l ixi,...,Xj^t, ...,Xn) = em,,N {xi, ...,Xn) - Xtem~l,N~l {xi, ■■■, Xj^.^, ..., Xn) ■ 



Wielding these two identities it is possible to build an 
efhcicnt algorithm to generate the full set of el for i ~ 
1, ...,N that is required to calculate the elements of the 
form Y^ ^{ui, Vi)Ji for the N x N Slater matrix. 



Once the Slater matrix is populated, the remainder 
of the work is involved in calculating its determinant. 
The full algorithm evaluates the probability density of 
a given CF state for a given set of co-ordinates. Based 
on the times recorded when running the program, our 



I 

full algorithm evaluation time scales roughly as iV^-^- 
A^^-^^, the better figure occurring for n = 1 and the worse 
figure for n = 4, which confirms our assertion that our 
algorithm is more efficient for smaller values of n. 

We note that the evaluation of the CF wavefunctions 
tends to suffer from a numerical precision issue, particu- 
larly for large N. To overcome the issue we simply store 
all numerical values to a higher precision — although this 
slows down the algorithm considerably, we are able to 
obtain accurate results up to around N = 40, which is 
high enough for our purposes. 



For classic reviews of quantum Hall physics, see R. Prange 
and S.M. Girvin, The Quantum Hall Effect (Springer- 
Verlag, New York, 1987). 

2 J. G. Groshaus et al, Phys. Rev. Lett. 98, 156803 (2007). 

^ I. V. Kukushkin, K. von Klitzing, and K. Eberl, Phys. Rev. 
Lett. 82, 3665 (1999). 

* P. Khandelwal et al, Phys. Rev. Lett. 81, 673 (1998). 

^ H. Cho et al., Phys. Rev. Lett. 81, 2522 (1998). 

® M. J. Manfra et al, Phys. Rev. B 54, 17327(R) (1996). 

'^ R. R. Du et al, Phys. Rev. Lett. 75, 3926 (1995). 

® J. K. Jain, Composite Ferm.ions, (Cambridge University 
Press, Cambridge, 2007). 

® K. Park and J. K. Jain, Solid State Commun. 119, 291 

(2001). 
"' J. K. Jain and R. K. Kamilla, Phys. Rev. B 55, R4895 

(1997). 
" G. MoUer and S. H. Simon, Phys. Rev. B 72, 045344 

(2005). 
^^ M. R. Peterson, T. Jolicoeur, and S. D. Sarnia, Phys. Rev. 

B 78, 155308 (2008). 
^■^ Note that the wavefunction of the form ip{ui,Vi) given 
in Eq. ((9]) is not quite complete: in fact wavefunctions 
$ {ui, Vi) describing A^ fermions must be fully antisymmet- 
ric under any exchange of particle labels, but ip {ui, Vi) in 
Eq. (O is only antisymmetric under exchange of labels 1 to 
A^^ and in labels A'^ -I- 1 to iVj,. The correct expression for 
the full wavefunction would be obtained by including a spin 
part to the wavefunction, X = |tt • • • 4^), and then fully 
antisymmetrizing by applying the antisymmetrizing oper- 
ator A i.e., we must construct ^{ui,Vi) = Aip {ui,Vi)X. 
In practice, when calculating matrix elements of the form 
J ^^0$dQ,i . . . dfijv, as long as the operator O does not 
depend on the spin state the matrix element reduces to 
NlJili'Otjj because the spin wavefunctions are orthonor- 



mal. Thus the version of the wavefunction given in Eq. ^ 
is sufficient. 

R. Morf and B. I. Halperin, Z. Phys B 68, 391 (1987). 
C. Toke, M. R. Peterson, G. S. Jeon, and J. K. Jain, Phys. 
Rev. B 72, 125315 (2005). 

F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 
The differences between the disc geometry case and the 
sphere geometry case are the values of the coefficients Ci in 
Eq. (|ll|l : in the disc geometry version of the potential the 
coefficients are derived using the disc geometry versions 
of the Haldane pseudopotentials, whereas in the sphere 
geometry the coefficients are derived using the sphere ge- 
ometry versions of the pseudopotentials (see, for example, 
Ref . [la ) . In the disc geometry the values of the coefficients 
are independent of the flux; however, this is no longer true 
in the sphere geometry and a different set of coefficients 
(and therefore a different effective potential) must be cal- 
culated for each value of the ffux. We have checked that in 
the thermodynamic limit the values of the sphere geome- 
try coefficients tend toward the same fixed values given in 
Table U] and thus it is sufficient to only consider the disc 
geometry version of the potential. 

C. Toke, and J. K. Jain, Phys. Rev. Lett. 96, 246805 
(2006). 

In principle, tilting the magnetic field can alter the effective 
interaction between the electrons by modifying the form of 
the wavefunction in the direction perpendicular to the elec- 
tron layer. We have not accounted for this effect. However, 
we note that the in-plane magnetic length is never much 
less than the well width so this effect is not expected to be 
substantial (see Table HV)) . 

F. Stern and W. E. Howard, Phys. Rev. 163 816 (1967). 
M. Stern et al, Phys. Rev. Lett! 108, 066810 (2012). 
L. Tiemann, G. Gamez, N. Kumada, and K. Muraki, Sci- 



18 



ence 335, 828-831 (2012). 

R. H. Morf, Phys. Rev. Lett. 80, 1505 (1998). 

E. H. Rezayi and S. H. Simon, Phys. Rev. Lett. 106, 

116801 (2011). 

See for example, M. E. J. Newman, and G. T. Barkema 

Monte Carlo Methods in Statistical Physics, (Clarendon 

Press, Oxford, 1999). 

K. Park, Ph.D. thesis. State University of New York at 



Stony Brook, 2000. 

L G. Macdonald, Symmetric Functions and Hall Polyno- 
mials, (Oxford University Press, Oxford, 1979). 
W.Pan, K.W. Baldwin, K.W. West, L.N. Pfeiffer, and D.C. 
Tsui, arXiv:cond-mat/1204.0557 (2012). 
R. Morf, and N. d'Ambrumenil, Phys. Rev. Lett. 74, 5116- 
5119 (1995). 



